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ABSTRACT 


The  California  l  ndercurrent  (CtJC).  a  poleward- flowing  feature  over  the  continental  slope,  is  a  key  transport 
pathway  along  the  west  coast  of  North  America  and  an  important  component  of  regional  upwelling  dynamics,  ‘litis 
study  examines  the  polewaid  undeicurrent  and  alongshore  pressure  gradients  in  the  northern  California  Current 
System  (ex's),  where  local  wind  stress  forcing  is  relatively  weak.  The  dynamics  oi  the  undercurrent  are  compared 
in  the  primitive  equation  Navy  Coastal  Ocean  Model  and  a  linear  coastal  trapped  wave  model.  Both  models  are 
validated  using  hydrographic  data  and  current-meter  oteervalions  in  the  core  of  the  undercurrent  in  the  northern 
ex's.  In  the  linear  model,  variability  in  the  predominantly  equatorwund  wind  stress  along  the  U.S.  West  Coast 
produces  episodic  reversals  to  poleward  How  over  the  northern  CCS  slope  during  summer.  I  lowever,  reproducing 
the  pc i -si si e nee  of  the  undercurrent  during  late  summer  requires  additional  incoming  energy  from  sea  level  var¬ 
iability  ajplied  south  of  the  region  of  the  strongest  wind  forcing.  The  relative  importance  of  the  burotropic  and 
bait'd inic  components  of  the  nicxJeled  alongshore  pressure  gradient  changes  with  latitude.  In  contrast  to  the 
southern  and  central  portions  of  the  OX'S,  the  barn  clinic  comjxment  of  the  alongshore  piessurc  gradient  provides 
the  primary  jxileward  force  at  Ct  C  depths  overthe  noilhent  (X^S  slope.  At  lime  scales  from  weeks  to  months,  the 
alongshore  pressure  gradient  force  is  primarily  balanced  by  the  Coriolis  force  associated  with  onshore  llow. 


1.  Introduction 


the  North  Pacific  Ocean,  transporting  heat,  chemical 
tracers,  and  organisms  over  alongshore  distances  of 
thousands  of  kilometers.  Relatively  warm,  saline,  and 
oxygen-depleted  water  is  associated  with  the  CUC 
from  Baja  California  to  Alaska  (Ilickcv  1979;  Thomson 
and  Krassovski  2010).  These  equatorial  water  properties 
are  transferred  into  the  interior  of  the  North  Pacific  by 
eddies  formed  from  instabilities  within  the  Cl  ?C  ( I  Iuyer 
et  al.  1998;  Garfield  <et  al.  1999).  The  CUC  modifies 
nutrient  concentrations  along  its  path  by  transporting 
water  with  low  ratios  of  nitrate  to  phosphate  (Liu  and 
Kaplan  1989;  Castro  ct  al.  2001)  and  also  influences  the 
distribution  of  marine  organisms  such  as  zooplankton 
(Swartz man  et  al.  2005).  Poleward  undercurrents  similar 


The  California  Undercurrent  (CUC)  flows  poleward 
over  the  continental  slope  along  the  eastern  boundary  of 
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to  the  Cl  'C  arc  common  in  eastern  boundary  current 
systems,  and  their  presence  can  greatly  affect  the  phys¬ 
ical  dy  namics  of  coastal  up  welling  (Hill  ct  al.  1  998), 

The  ('l  C  is  a  seasonal  feature  in  the  northern  California 
eurrent  System  (CX'S:  Hickey  l1)?1):  Thomson  and 
Krassovski  2010).  which  includes  the  coastal  waters 
offshore  of  northern  California  (CA).  Oregon  (OR). 
Washington  (WA).  and  southern  British  Columbia 
(BC).  Following  Hickey  (1979).  the  CliC  is  de lined 
as  a  subsurface  maximum  in  poleward  How  over  the 
continental  slope  and  below  the  main  pycnocline.  At 
these  latitudes,  the  CCChasa  typical  volume  transport 
of  0.5-1.5Sverdrups  (Sv  1  S\  =  I06m3s  l).  with  a 
(U-20-km-wide  subsurface  core  over  the  upper  slope 
(Ilickex  1979;  Pierce  et  al  2(KK)).  In  the  seasonal  cycle, 
poleward  How  appears  over  the  upper  slope  at  depths 
<3(X)m  during  July  and  strengthens  through  late  October 
(I  lickex  1979;  Thomson  and  Krassovski  2010).  Digrangian 
floats  at  250-600-m  depth  between  40  and  47  N  indicate 
that  poleward  How  is  present  year-round  but  is  weakest 
during  March  and  strongest  during  October  (Collins  et  al 
2003).  As  the  up  welling  season  progresses  from  late  April 
to  early  October,  the  core  of  the  undercurrent  shoals 
from  below  500-m  depth  to  the  upper  slope  (Pelland  et  al. 
2013).  and  the  relatively  warm  and  saline  Ct‘C  water 
mass  forms  a  greater  percentage  of  water  over  the  shelf 
(MacFadyen  et  al.  2(X)S)  During  winter,  winds  are  down- 
welling  favorable,  and  the  CliC  merges  with  the  surface- 
in  tens  i  tied  poleward  How  over  the  slope  commonly  known 
as  the  Davidson  Current  (Hickey  1979).  Ilie  spring  tran¬ 
sition  to  equatorward  How  over  the  upper  slope  precedes 
the  transition  t  o  upxx  el  ling-fax  arable  wind  stress  by  around 
one  month  ( ITiomson  and  Krassovski  2010). 

Ilie  CliC'  influences  the  seasonal  dynamics  of  coastal 
upwclling  in  the  CCS  because  of  its  relationship  to  the 
alongshore  pressure  gradient  and  the  bottom  boundary 
layer.  In  the  northern  CCS.  the  baro tropic  component 
of  the  alongshore  pressure  gradient  force  (API  )  at  the 
coast  is  directed  equatorward  during  winter  and  spring 
and  becomes  poleward  during  the  summer  upwclling 
season  (Hickey  and  Pola  1983).  Imposing  a  poleward  API 
in  two-dimensional  models  of  coastal  upwclling  produces 
a  poleward  undercurrent  over  the  slope  and  shifts  the 
depth  of  onshore  How  from  the  bottom  boundary  layer 
to  the  interior  (Werner  and  Hickey  19, S3;  Federiuk  and 
Allen  1995).  Similarly.  Pringle  and  Dever  (2009)  found 
that  incorporating  a  realistic  CliC  into  a  three-dimensional 
regional  hindeast  model  leads  to  shallower  source  depths 
for  upwclling.  Idealized  modeling  studies  of  upwclling 
ecosystems  and  biogcochemistiy  often  use  an  externally 
applied  poleward  API  to  incorporate  undercurrent  dy¬ 
namics  (Laihuiliere  et  al.  2010;  Siedlecki  et  al.  2012). 
I  lowever.  it  is  not  clear  how  the  structure  of  the  API  varies 


in  the  cross-shore  direction  or  whether  there  is  a  baroelinic 
component  associated  with  large-scale  density  gradients. 

Despite  the  prevalence  and  importance  of  eastern 
boundary  poleward  undercurrents  like  the  CIX\  mech¬ 
anisms  for  their  generation  are  not  fully  understood. 
Hill  et  al.  (1998)  classify  three  proposed  mechanisms: 
1 )  a  response  to  wind  stress  near  the  coast.  2)  pressure 
gradients  originating  in  the  open  ocean,  and  3)  recti ti- 
cation  of  oscillating  flow  over  complex  topography .  The 
first  two  mechanisms  depend  on  a  large-scale  poleward 
API*.  In  idealized  wind-driven  models  with  stratifica¬ 
tion.  a  poleward  API  and  undercurrent  can  be  gener¬ 
ated  if  the  upwclling-fav orable  winds  hav  e  an  alongshore 
structure,  an  idea  lirst  recognized  by  Yoshida  (1967). 
However,  an  open-ocean  pressure  gradient  associated 
with  a  poleward  decrease  in  temperature  and  sea  surface 
height  (SSI  I )  can  also  be  associated  with  a  poleward  flow 
oxer  the  slope  (Iiuthnancc  1984:  Csanndv  1985).  In  the 
third  mechanism,  polcw  ard  How  is  associated  with  a  mean 
asymmetrical  distribution  of  pressure  across  topographic 
features  (Holloway  1987:  Brink  2010.201 1 ).  In  this  case, 
the  oscillating  flow  may  be  forced  by  variable  alongshore 
wind  stress,  but  oxer  the  course  of  many  oscillations, 
mean  poleward  flow  oxer  the  slope  arises  in  the  absence 
of  a  mean  wind  stress.  All  of  the  above  mechanisms  can 
produce  poleward  How  oxer  the  slope  that  qualitatively 
resembles  the  cross-shelf  and  depth  structure  found  in  the 
linute  d  existing  observations  of  the  undercurrent,  but  the 
different  types  of  forcing  should  result  in  differences  in 
the  temporal  variability  of  poleward  How.  The  absence  of 
a  clearly  defined  mechanism  for  undercurrent  generation 
continues  due  to  the  low  number  of  velocity  time  series 
characterizing  the  temporal  variability  of  slope  currents 
and  the  difficulty  in  obtaining  measurements  of  the  APT. 

In  the  theoretical  work  that  emphasizes  wind  forcing 
as  a  mechanism  for  undercurrent  generation,  coastal 
trapped  waxes  (CTWs)  play  an  important  role.  (TWs 
arc  hybrids  of  internal  Kelvin  waxes  and  harotropic 
continental  shelf  waxes  (Wang  and  Mooers  1976).  In 
pure  internal  Kclxin  waves*  over  a  flat  bottom,  along¬ 
shore  pressure  gradients  accelerate  alongshore  motions 
near  a  coastal  wall  but  cross -shore  velocity  remains  zero. 
In  shelf  waxes,  cross- shore  motions  act  to  conserve  po¬ 
tential  vorticity  oxer  a  sloping  bottom.  When  stratifica¬ 
tion  and  bottom  slope  are  present.  (TWs  arc  nearly 
harotropic  oxer  the  shelf,  xxith  a  more  barOclinic  struc¬ 
ture  over  the  slope.  Higher  modes  tend  to  have  more 
baroelinic  structure  than  lower  modes. are  damped  more 
rapidly  hx  friction,  and  propagate  more  slowly. 

I'he  undercurrent  is  a  robust  feature  of  models  forced 
by  wind  stress  in  a  region  of  limited  alongshore  extent. 
In  stratified  models  with  a  flat  bottom.  Kelvin  waves  are 
generated  at  the  southern  edge  of  the  forcing  region,  and 
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a  poleward  APF  that  accelerates  the  undercurrent  is  set 
up  in  their  wake  (McCreary  llJSI;  Philander  and  Yoon 
PJS2).  When  a  continental  shelf  is  present,  the  first- 
mode  (TW  is  associated  with  cquatorward  flow,  and 
a  poleward  flow  develops  over  the  slope  with  the  passage 
of  a  second-mode  (1  W  (Suginohara  ll)74,  P)S2). 

Hie  two  essential  ingredients  for  wind -driven  poleward 
undercurrents  are  alongshore  structure  in  the  wind  forc¬ 
ing,  which  generates  an  APF.  and  stratification,  which 
allows  vertical  shear  in  the  alongshore  currents.  With 
these  basic  factors  included,  the  two-layer  analytical 
model  of  Yoshida  (P)SO)  produces  a  poleward  under¬ 
current  over  nearly  the  entire  range  of  length  scales  and 
frequencies.  ITie  b  effect  is  not  essential  for  modeling  the 
undercurrent,  but  offshore  Rossbv  wave  propagation  can 
lead  to  a  shallowing  and  intensification  of  the  undercurrent 
(Suginohara  and  Kitamura  P1S4:  McCreary  and  Chao 
P)S5:  Marchesiello  et  al.  2003).  Positive  wind  stress  curl 
near  the  coast  on  a  b  plane  can  also  strengthen  the  un¬ 
dercurrent  (McCreary  and  Chao  FKS;  Batteen  1W7). 

Ilie  timing  of  Cl;Cde\  elopment  in  the  northern  CCS. 
shown  by  Thomson  and  Krassovski  (201b).  is  qualita¬ 
tively  consistent  with  a  summertime  poleward  APF  in¬ 
ferred  from  coastal  sea  level  observations  (Hickey  and 
Pola  F)S3).  I  Iickcy  and  Pola  (P)S3)  attributed  the  sea 
level  gradient  in  the  northern  CCS  to  the  alongshore 
structure  of  wind  stress  using  the  steady  baro tropic 
model  of  Csanady  (  N7S)/I  his  structure  in  summer  wind 
stress,  with  a  maximum  off  of  northern  CA.  has  been 
linked  to  the  Cl’C  and  poleward  APF  in  the  baroclinic 
models  (e.g.,  McCreary  P)S1)  described  above.  How¬ 
ever.  direct  comparisons  with  observed  variability  in  the 
(VC  have  been  limited,  making  it  difficult  to  rule  out 
other  mechanisms  of  producing  poleward  llow.  Fur¬ 
thermore.  because  the  API  is  difficult  to  measure  away 
from  the  coast,  its  structure  and  variability  over  the 
slope  remain  poorly  understood. 

The  goals  of  this  paper  are  to  investigate  the  spatial 
and  temporal  structure  of  the  poleward  undercurrent 
and  the  alongshore  pressure  gradient  along  the  U.S. 
West  Coast,  extending  previous  work  by  incorporating 
realistic  wind  forcing  and  boundary  conditions  at  time 
scales  from  weather  events  to  seasons.  Emphasis  is 
placed  on  the  development  of  the  CUC  in  the  northern 
portion  of  the  CCS.  where  local  wind  stress  forcing  is 
relatively  weak.  The  analysis  uses  two  numerical  hind- 
casts  with  different  levels  of  complexity.  The  Navy 
Coastal  Ocean  Model  of  the  CCS  (NCOM-CCS)  is  used 
to  more  fully  understand  the  APF  over  the  continental 
slope,  where  it  is  difficult  to  observe.  A  linear  Cl  W 
model  is  also  used,  which  is  forced  by  temporally  and 
spatially  variable  winds  along  the  coast  hut  excludes  the 
effects  of  complex  topography  ,  open-ocean  forcing,  wind 
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stress  curl,  and  the  b  effect.  These  limitations  are  less 
restrictive  in  the  northern  CCS  than  at  locations  farther 
south  because  the  region  lies  north  of  latitudes  where 
seasonal  Rossbv  wave  propagation  occurs  (Clarke  and 
Shi  llWl:  Kelly  et  al.  PW3)  and  wind  stress  curl  is  rela¬ 
tively  weak  (Bakun  and  Nelson  llMl).  I  his  model  can  be 
directly  compared  with  observations  and  includes  the 
basic  ingredients  in  wind-driven  theories  of  undercurrent 
generation:  alongshore  structure  in  wind  stress,  along¬ 
shore  pressure  gradients,  and  stratification.  Ilie  more 
complete  physics  of  NCOM-CCS  provide  our  best  esti¬ 
mate  of  the  ocean  state  at  locations  where  observations 
are  limited,  while  the  CTW  model  is  used  to  isolate  the 
linear  response  to  alongshore  wind  stress  in  the  CCS. 

Both  models  are  compared  with  an  extensive  array  of 
observations,  including  currents  from  moored  sensors  and 
hydrographic  data  (section  3).  ITie  (TW  model  is 
adjusted  to  show  the  scnsiti\  itv  of  the  alongshore  velocity 
to  friction,  number  of  modes.  local  wind  stress,  and  the 
southern  boundary  condition  (section  4a).  In  this  linear 
model,  alongshore  winds  along  the  U.S.  West  Coast  gen¬ 
erate  episodic  periods  of  poleward  flow  over  the  northern 
CCS  slope  during  summer,  when  wind-driven  flow  over 
the  shelf  is  predominantly  cquatorward.  However,  the 
inclusion  of  sea  level  data  at  the  southern  boundary. 
>13(K)km  to  the  south,  is  required  to  reproduce  the  ob¬ 
served  strength  and  persistence  of  the  slope  undercurrent 
during  late  summer.  The  APF  is  examined  near  the  coast 
and  over  the  slope  in  both  models  (section  lb),  then 
connected  to  the  dynamics  of  Cl  IC  development  through 
the  analysis  of  the  alongshore  momentum  balance 
(section  4c).  It  is  shown  that  alongshore  density  gradients 
contribute  to  a  poleward  API  over  the  northern  CCS 
slope  in  contrast  to  locations  farther  south. 

2.  Methods 

a.  Observations 
l)  Coastal  sea  level 

Coastal  sea  le\el  time  series  were  obtained  for  M)  tide 
gauge  locations  (Fig.  Ia).  Sea  level  data  were  adjusted 
by  adding  equivalent  sea  lex  cl  pressure  (0.01  m  mbar  1 
atmospheric  pressure)  measured  at  nearby  stations 
(Fig.  la)  in  order  to  estimate  subsurface  pressure.  I  Iourly 
time  series  w  ere  low-pass  filtered  using  a  cosine  Lane/ os 
filter  with  a  half-power  point  of  46  h  and  decimated  every 
6h.  Because  the  absolute  vertical  displacement  from 
a  geopotential  surface  is  unknown,  mean  values  are  re¬ 
moved  so  that  time  scries  represent  anomalies.  Tide 
gauge  data  were  obtained  from  the  U.S.  National  Ocean 
Service,  and  atmospheric  pressure  data  were  obtained 
from  the  National  Buoy  Data  Center. 
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NCOM  SSH  (m)  NCOM  depth  avg  v  (m/s)  CCMP  10m  wind  (m/s) 


128  126  124  122  120  118  116  128  126  124  122  120  118  116  128  126  124  122  120  118  116 


Fir,.  I.  Mean  model  fields  (NCOM-CC'S.  with  global  NCOM  north  of  49  N)  and  atmospheric  forcing  (CC’MP)  for  the  California 
Current.  15  Jun  15  Sep  2005.  The  200- .  500- .  and  1000-m  isobaths  are  shown  as  black  contours,  (a)  SSH.  with  tM)5-ni  coni  our  intervals  in 
while.  Red  squares  indicale  tide  gauges  ai  San  Diego  (SID).  Port  San  Luis  (PSL).  Monterey  I3av  (MB).  Point  Reyes  (PR).  Arena  Cove 
(AC).  Humboldt  Ray  (HB),  Newport  (Nil).  La  Push  (LP).  and  Neah  Bay  (NB).  Red  circles  show  nearby  measurements  of  sea  level 
pressure,  (b)  Northward  component  of  depth -averaged  velocity.  Gray  triangles  indicate  mooring  locations,  (c)  Wind  speed  and  velocity 
vector.  Gray  lines  indicate  cross-shore  segments  used  in  the  coastal  trapped  wave  model.  (TW  model  segments  coincide  with  tide  gauge 
locations  at  the  coast,  except  for  the  WA  and  To  lino  (TF)  segments.  Geographic  boundaries  of  BC,  WA.  OR.  and  ('A  are  indicated. 


2)  Mooring  time  series 

Velocity  data  arc  u^d  from  two  locations,  one  on  the 
British  Columbia  slope  and  one  on  the  Washington  shelf 
(Fig  lb).  Hie  Al  mooring,  where  data  have  been  col¬ 
lected  from  1 0S5  through  the  present,  is  located  on  the 
>00-m  isobath  at  4K°32f\,  126  12AV  flTtomson  and 
Krassovski  2010).  During  2001  and  2(X).\  velocity  data 
were  collected  at  nominal  depths  of  35.  1(X ).  175,  and 
3(X)  m  using  Aanderaa  recording  current  meters  (RCM)*S, 
although  actual  current- me  ter  depth  and  data  availability 
vary  between  deployments.  Velocity  components  were 
rotated  25  counterclockwise*  parallel  to  the  local  isobath 
orientation,  and  time  series  were  filtered  using  a  Kaiser- 
Bessel  low-pass  filter  with  a  30-h  cutoff  and  decimated  to 
daily  intervals.  On  the  Washington  shelf,  maintained  as 
part  of  the  River  Influences  on  Shelf  Ecosystems  ( RISK) 
project  (Mickey  et  al.  2010).  the  RISK  North  (RN) 
mooring  is  located  on  the  70- m  isobath  at  47  (X)-N. 
124  30-W  (data  provided  by  H.  Dcvcr  2(X)6,  personal 
communication).  Velocity  components  at  the  RN  site 


were  rotated  10  clockwise*  and  time  series  were  filtered 
and  decimated  in  the  same  manner  as  coastal  sea  level. 

3)  GEOSTROPI IK’  VELOCITY 

Conductivity',  temperature,  and  depth  (CTD)  mea¬ 
surements  were  obtained  during  September  2CX)5  as  part 
of  the  Ecology  and  Oceanography  of  Harmful  Algal 
Blooms  Pacific  Northwest  (ECOHAB-PNW)  project 
(MacFadven  et  al.  2C XKS).  Gcostrophic  velocity,  refer¬ 
enced  to  500dbar.  was  calculated  using  the  method  of 
Reid  and  Mantyla  ( F)76),  Because  shipboard  ADCPdata 
are  not  available  at  5(X)dbar  to  compute  absolute  ve¬ 
locity.  NCOM-CCS  velocity  at  5(H) dbar  is  used  as  a 
reference  for  consistent  comparison  between  model 
and  observations. 

h.  Mfulcls 

1 )  Navy  c  oastal  oc  ean  model 

NCOM-CCS  (Shulman  et  al.  2C X )7 )  is  a  primitive  equa¬ 
tion.  3D,  hydrostatic  model  with  a  horizontal  resolution  of 
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-9  km.  I  hc  model  domain  covers  the  region  30  — 19  N. 
135^-115  W  (Figs.  la.b).  and  hindcast  results  are  pre¬ 
sented  for  the  year  2005.  The  model  is  forced  with 
atmospheric  products  derived  from  the  Coupled  Ocean- 
Atmosphere  Mesoscale  Prediction  System  (COAMPS: 
Iloelur  1997).  I  he  model  assimilates  satellite -derived 
SSI  I  and  sea  surface  temperature  (SST)  data  via  the 
Modular  Ocean  Data  Assimilation  System  (MOD AS: 
Fox  et  al.  2002).  which  uses  the  SST  and  SSI  I  data  to 
generate  synthetic  temperature  and  salinity  profiles. 
Open  boundary  conditions  for  the  regional  N  COM -CCS 
arc  derived  from  the  global  NCOM  (Rhodes  et  al.  2(X)2; 
Barron  et  al.  2004).  which  has  a  Vh  horizontal  resolution. 
In  the  version  used  here,  the  vertical  coordinate  system 
is  composed  of  40  levels  total.  19  s  (terrain  following) 
levels  on  top  of  21  c  (constant  depth)  levels  (Shulman 
et  al.  2(K >7).  There  is  no  tidal  forcing  or  river  input,  and 
the  mouth  of  the  Strait  of  .luan  de  Fuca  at  Noth  Bay 
(NB:  I'ig.  la)  is  closed  off  by  a  straight  coastline.  I  se  of 
c  levels  below  138  m.  which  is  typically  shallower  than 
the  depth  of  the  shelf  break  in  the  CCS.  avoids  errors  in 
calculation  of  the  pressure  gradient  that  can  arise  near 
the  steep  topograph}  of  the  continental  slope  when  us¬ 
ing  s  coordinates  (Ilanev  1991). 

2)  Coastal  trapped  wave  model 

A  (TAV  model  is  used  to  further  isolate  mechanisms 
for  subinertial  variability  in  velocity  and  sea  level  along 
the  west  coast  of  North  America.  Derivation  and  appli¬ 
cation  of  CTVV  theory  are  discussed  extensively  by  Clarke 
and  Van  Ciorder  (1986)  and  will  only  be  summarized 
briefly  here.  'Time-dependent  effects  of  alongshore  wind 
stress  are  included  in  the  model,  along  with  rotation, 
continuous  stratification*  bottom  slope,  and  weak  bottom 
friction.  Wind  stress  curl  and  the  b  effect  are  neglected. 
'Hie  flow  is  assumed  to  be  linear  and  Boussinesq.  with 
gradually  v  a  lying  topography  in  the  alongshore  direction. 

The  coordinate  system  is  aligned  such  that  a  is  the 
cross-shore  distance  from  the  coast  (positive  onshore), 
y  is  the  alongshore  position  (positive  northward),  and 
z  is  the  vertical  position  (positive  upward).  The 
problem  is  greatly  simplified  in  the  long-wave  limit,  in 
which  cross-shelf  scales  arc  assumed  to  be  much 
shorter  than  alongshore  scales  (i.e.,  *7  y"  "7  a~). 
and  time  scales  of  variability  are  assumed  to  be  much 
longer  than  the  inertial  period  (i.e.,  *7  r  /".where 
/  is  the  Coriolis  frequency).  Turbulent  stresses  are 
restricted  to  infinitesimally  thm  boundary  layers  at  the 
surface  and  bottom.  In  terms  of  pressure,  the  system 
can  be  written  as 

PM,  +  (.1) 


where  A*(r)  is  the  buoyancy  frequency  of  the  back¬ 
ground  state,  and  p  is  the  deviation  of  pressure  from  the 
background  state.  Boundary  conditions  are  give n  by 


A’2 

pm  — p  —  0  at  r-0. 

(2) 

rp.  .  ft' 

,,  at  ,-0. 

(3) 

I‘X  (/\ ,  f  fPy  )  +  ( t‘X  ’  P  X  -  -  (  1 

at 

//(a),  and 

(4) 

px—  0  as  a/ 

(5) 

where  t*v(y,  t)  is  the  alongshore  component  of  wind 
stress,  //(a*.  \ )  is  bottom  depth  (assumed  to  vary  gradu¬ 
ally  in  the  alongshore  direction),  r  is  a  linear  friction 
coefficient  such  that  the  bottom  stress  t£  -r0/*y.  r0  is 
a  reference  density,  and  y  is  alongshore  velocity.  Ill  esc 
boundary  conditions  represent  a  free  surface  |(2)|.  no 
net  cross-shelf  transport  at  the  coastal  wall  |(3)|.  Iikman 
transport  governed  by  no  normal  flow  at  the  bottom 
[(4)|.  and  coastal  trapping  |(5)|.  Solutions  for  pressure 
arc  separated  into  orthogonal  components. 


/>(.v.  v.r.f)  -  □  f'nix.  r)f  ;,(y.f). 


((>) 


n  i 


where  I:n  is  the  frictionless  free-wave  structure  of  the  nth 
mode,  and  f  „  is  the  amplitude  of  the  //th-mode  response. 
The  f  /r  satisfy 


□ 


□  t— * 1  D  (itVH  '  #W 
m  I 


(7) 


where  cn  is  the  phase  speed  of  the  nth  mode.  hn  is  a  wind 
coupling  coefficient.  a,w  is  a  frictional  decay  cocfticient. 
and  aMn are  coefficients  for  frictional  coupling  to  the  mth 
mode.  Numerical  methods  forcalculating  the  five-wave 
properties  f:„.  c„.  hn.  and  a/w,  from  mean  stratification 
A:2(r).  cross-shore  bottom  topography //(.v).  and  Coriolis 
parameter  f  are  described  by  Brink  and  Chapman 
(1987).  Versions  of  the  software  for  the  Matrix  labo¬ 
ratory  (MAT LAB)  programming  environment  were 
provided  bv  Dr.  Kenneth  Brink  [Woods  Hole  Oceano¬ 
graphic  Institution  (WI IOI )].  Wind  coupling  coefficients 
hn  and  frictional  coefficients  a/v1,  are  calculated  using  the 
energy -conserving  normalization  of  Brink  (1989).  Once 
the  free-wave  properties  are  found,  solutions  to  the 
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coupled  equations  in  (7)  can  be  obtained  using  the  nu¬ 
merical  method  of  characteristics  described  by  Oarkc 
and  Van  Go rde r  ( 1 0 86 ) . 

The  free -wave  parameters  were  calculated  for  10 
coastal  segments  (Fig.  1c).  Mean  A"(r)  profiles  for  each 
segment  were  computed  from  C I  D  data  in  the  World 
Ocean  Database  (Boyer  et  al.  2(X)6).  using  casts  offshore 
of  the  KMX  km  isobath  but  no  farther  than  40()-km  offshore. 
Previous  studies  have  shown  that  varying  A;'(r)  produces 
small  (<10%)  changes  in  /■„,  c„.  and  hn  in  the  long-wave 
limit  but  can  change  u,u„  by  up  to  l(X)%  (Battisti  and 
1  tic  key  1984:  (Tiapman  1987).  Bottom  topography  is  from 
the  National  Geophysical  Data  (enter  coastal  relief 
mod  el  ( http://www.ngdc.  noaa  .gov/mgg/coastal/cr m.html ). 
Wind  stress  was  calculated  using  the  formula  of  Dirge  and 
Pond  ( l  DS l )  from  wind  velocity  interpolated  to  the  in¬ 
shore  50km  of  each  segment  from  the  Cross-Calib rated 
Multi-Platform  {('CMP)  product  (Fig.  Ic),  which  incor¬ 
porates  data  from  the  Quick  Scatterometer  (QuikSCAT) 
and  coastal  buoys  (Atlas  et  al.  201 1). 

To  resolve  the  slowest  phase  speed  in  the  model 
(0.39  ni  s  1  for  the  fourth  mode  off  of  southe rn  California ). 
Ay  and  Af  were  set  to  2  5  km  and  6h.  respectively.  ITie 
(TW  model  was  run  for  the  years  2(X)3-05.  with  results 
from  2(X)3  discarded  as  spinup.  Unless otherw ist.  stated, 
the  infinite  sums  in  (6)  and  (7)  were  approximated  using 
four  modes.  A  constant  linear  friction  coefficient  r  — 
2.5  x  10  4  ms  1  was  used  for  the  calculation  of  a,W}. 
ITi is  is  equivalent  to  the  value  used  for  the  continental 
slope  but  less  than  that  used  for  the  shelf  by  Chapman 
(1987).  At  the  southern  boundary,  f  i(0.  t)  -  r0gh/ 
f' |((k  0),  and  adjusted  sea  le\el  at  San  Diego  was  used 
for  h(y).  This  relationship  assumes  that  the  San  Diego 
slu  level  is  dominated  by  the  first  mode.  Sensitivity  to 
the  number  of  modes,  friction  coefficient,  presence  of 
local  wind  stress,  and  the  use  of  sea  level  data  at  the 
southern  boundary  will  be  explored. 

3.  Model  validation 

Validation  wi  th  observ  ations  is  necessary  before  using 
the  models  to  extend  prior  work  on  the  CFC  dynamics. 
The  models  are  first  compared  with  observed  time  scries 
of  alongshore  velocity  at  the  RN  mooring  over  the  WA 
midshelf  (Fig.  lb),  where  previous  studies  have  dem¬ 
onstrated  the  importance  of  CTW  dynamics  (Battisti 
and  Hickey  1984).  To  assess  the  variability  and  strength 
of  the  modeled  CUC.  the  models  are  also  compared  w  ith 
time  series  at  the  Al  mooring  over  the  British  Columbia 
slope  (Fig,  lb). 

Statistical  agreement  between  the  models  and  obser¬ 
vations  is  quantified  using  both  correlation  coefficients 
(CC)  and  Willmott  skill  (WS).  Correlation  provides  a 


test  for  statistically  significant  agreement  at  weather- 
band  time  scales.  CC  were  calculated  from  a  detrended 
time  series,  and  significance  levels  were  determined  us¬ 
ing  effective  degrees  of  freedom  { Finery  and  Thomson 
2004.  p.  260).  WS  normalizes  the  mean  squared  error, 
which  includes  bias  and  differences  in  variance  and 
trends,  allowing  for  comparisons  of  skill  between  dif¬ 
ferent  models  at  multiple  observational  locations.  WS 
v  aries  between  0  and  1  and  is  defined  as 


~  2  *  V"/ 

( in  o  o  o  ) 

whe  re  ///  is  a  modeled  \  ariablc.  o  is  an  observed  \  ariable. 
and  the  angle  brackets  denote  an  average  over  the  series 
(Willmott  1981). 

a.  Vitlshclf 

During  June-September  2005.  observed  currents 
were  predominantly  cquatorward  in  the  middle  of  the 
water  column  over  the  WA  shelf  (Fig.  2a),  and  sig¬ 
nificant  variability  was  also  observed  at  weather-band 
( -5  days)  time  scales  (Fig.  2a).  Reversals  to  poleward 
flow  do  occur  during  this  period  over  the  WA  shelf, 
becoming  more  frequent  during  September.  Correla¬ 
tions  between  NCOM-CCS  and  observations  at  RN  are 
similar  throughout  the  water  column  (CC  -  0.68-0.74), 
while  CTW  model  correlation  is  strongest  in  the  lower 
half  of  the  water  column  (CC  —  0.61-0.66)  (Fig.  3a). 
All  correlations  arc  significant  at  the  95%  confidence 
level.  A  similar  increase  with  depth  in  the  correlation 
between  a  Cl  W  model  and  obse  mil  ions  was  reported 
by  Chapman  (1987)  over  the  California  she  If.  Profiles  of 
WS  from  the  two  different  models  indicate  that  the  Cl  W 
model  generally  has  weaker  skill,  particularly  in  the  lower 
part  of  the  water  column  (Fig.  3a).  TTic  CTW  model  most 
accurately  predicts  the  timing  of  fluctuations  near  the 
bottom,  but  the  magnitude  of  the  modeled  velocity 
agrees  more  closely  with  observations  near  the  surface. 

h.  Slope 

During  2(X)5  at  the  Al  site  on  the  continental  slope, 
models  show  greater  agreement  with  observations  at 
undercurrent  depths  than  near  the  surface  (Figs.  2b.c). 
Observations  at  35 m  show  cquatorward  flow  that  in¬ 
tensifies  from  June  to  September  (Fig.  2b),  a  pattern  that 
is  not  shown  in  either  model.  At  300-m  depth,  the 
strongest  cquatorward  flows  were  observed  during  late 
May.  following  the  transition  in  local  wind  stress  (Fig.  2c). 
By  mid-Julv.  the  flow  is  predominantly  poleward  but 
with  considerable  event-scale  variability  throughout  the 
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NCOM  and  CTW  models.  RN  mooring  velocity  35  m  {70-m  isobath)  -  2005 


NCOM  and  CTW  models.  A1  mooring  velocity  35  m  {500-m  isobath)  -  2005 
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NCOM  and  CTW  models.  At  moonng  velocity  300  m  (500-m  isobath)  -  2005 
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Fill.  2.  Comparison  of  CTW  and  NCOM-CCS  modeled  alongshore  velocity  with  cuirent- 
meier  measurements  at  (a)  35-m  depth  at  the  WA  shelf  RN  mooring  siie  over  the  70-m  isobath, 
(b)  35-m  depth  at  the  I3C  slope  A 1  mooring  site  over  tltL  500-m  isobath,  and  (c)  300- m  depth  at 
the  I3C  slope  A  I  mooring  site  over  the  500-m  isobath.  Note  lhai  a  different  scale  has  been  used 
in  (c)  to  show  the  smaller  fluctuations  in  velocity  at  that  location. 


time  scries.  Both  models  show  a  transition  from  equa¬ 
tor-  to  poleward  flow  at  this  depth  during  summer  2005. 
although  the  (T  W  model  tends  to  overestimate  the 
magnitude  of  the  equatorward  flow.  During  2004,  the 
CTW  time  series  indicate  a  similar  seasonal  progression 
as  in  2005.  with  the  strongest  equatorward  flow  observed 
during  April  at  l(X)m  (Fig.  4).  However,  at  35m.  ob¬ 
served  equatorward  velocities  exceed  0.2  ms  1  during 
July,  which  is  not  reproduced  in  the  (TW  model. 

Over  the  5(Xkm  isobath,  the  largest  correlations  occur 
for  both  models  at  l(X)m  (Fig.  3b).  All  arc  significant  at  the 
95%  confidence  level.  Hie  CC  and  WS  statistics  both  in¬ 
dicate  that  NCOM-CCS  outperforms  the  CTW  model 
at  300  m.  If  comparisons  arc  restricted  to  the  same  time 
period  as  the  2005  RN  shelf  measurements  discussed  in 
the  previous  section,  the  correlation  with  observations 


improves  for  both  models  at  all  depths,  particularly  for 
NCOM-CCS  (Fig.  3c).  If  time  series  are  bandpass  filtered, 
using  a  Hanning  window  to  remove  periods  >20  days, 
correlation  coefficients  improve  substantially  for  the  CTW 
model  over  the  slope  (0.64. 0.79.  and  0.60  for  35.  1(XL  and 
300m.  respectively:  not  shown),  indicating  strong  agree¬ 
ment  at  weather-band  time  scales.  'I*hc  WS  metric  indicates 
that  both  models  have  relatively  strong  skill  ( —US)  at  300  m 
and  weak  skill  (  -0.3)  at  35  m.  The  weak  skill  at  35m  in¬ 
dicates  that  the  models  do  not  accurately  represent  the 
strength  and/or  cross-shone  position  of  the  equatorward 
shelf  break  jet  at  the  A l  site.  At  300m.  the  modeled  sea¬ 
sonal  progression  from  equator-  to  poleward  flow  results  in 
higher  skill  than  at  the  surface  during  this  time  period. 

A  hydrographic  section  offshore  of  Copalis  Beach. 
WA  ( 47HS-N)  shows  the  spatial  distribution  of  alongshore 


depth  (m)  depth  (m) 
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RN  mooring  (70~m  isobath),  21  May  -  05  Oct  2005 


A1  mooring,  21  May  -  05  Oct  2005 


FlO,  3,  Profiles  of  model  validation  statistics  for  (a)  the  RN 
mooring  over  the  70-m  isobath,  (h)  the  A I  mooting  over  the  500-nt 
isobath  for  2004  and  2005.  and  (c)  the  A I  mooting  restricted  to 
21  May  5  Oct  2005  (the  deployment  period  of  the  RN  mooting). 
Black  lines  are  CC;  gray  lines  are  WS  scores, 

velocity  on  13  September  2005  (Fig.  5a),  Over  the 
slope,  the  observed  core  of  the  Cl  IC  is  located  at  200  m. 
After  referencing  the  observ ations  to  0.05  m  s  1  at  5(X)dbar 
for  consistency  with  NCOM-CCS  (see  section  3).  the 


velocity  in  the  undercurrent  core  is  found  to  be 
(UWms  Equatorward  flow  is  present  over  the  outer 
shelf  and  slope  and  is  strongest  near  the  surface  over 
the  100-m  isobath.  Although  poleward  How  is  present 
over  the  slope  in  both  models  (figs.  5b.c).onlv  NCOM- 
CCS  shows  the  appropriate  distribution  of  equator- 
ward  flow.  'Fhe  CTW  model  shows  no  evidence  of  the 
"“0.3ms  1  near-surface  jet  that  extends  over  the  WA 
outer  shelf  and  slope,  unlike  NCOM-CCS  (Figs.  5h.c). 
I  he  location  of  maximum  equatorward  How  at  the 
surface  has  been  observed  to  migrate  offshore  from 
spring  to  summer  (Hickey  1989).  a  process  that  likely 
cannot  be  captured  by  the  CTW  model  that  neglects 
the  horizontal  advection  of  density  and  momentum. 
Weak  skill  at  35  m  over  the  British  Columbia  slope  in 
both  models  (I  igs,  2b.  3c.  and  4a)  may  be  related  to  the 
complex  topography  of  La  Perouse  Bank  located  on¬ 
shore  of  the  mooring  site,  as  well  as  CTW  scattering 
where  the  shelf  widens  from  WA  to  BC  (Wilkin  and 
Chapman  1990:  Hickey  et  al.  PHI).  'Hie  CTW  model 
does  not  include  scattering  by  alongshore  variations  in 
topography,  and  this  process  may  not  be  accurately 
resolved  bv  the  9-km  resolution  of  NCOM-CCS, 

Because  the  alongshore  currents  vary  on  the  time 
scales  of  weather  events,  synoptic  sections  of  geo- 
strophic  velocity  do  not  necessarily  represent  the  sea¬ 
sonal  mean.  Direct  current-meter  measurements  across 
the  Washington  shelf  and  slope  during  the  period 
23  Julv-26  August  1972  (Hickey  1979)  provide  an 
opportunity  for  comparison  of  the  model  results  with 
observations  of  the  mean  latc-summer  structure  of 
alongshore  velocity  over  the  shelf  and  slope,  although 
the  obse nations  arc  from  a  different  year.  The  un¬ 
dercurrent  is  clearly  present  in  NCOM-CCS  results  from 
the  same  part  of  the  season  (Fig,  5d).  but  the  0.09ms  1 
CUC  core  at  500-400  m  is  both  deeper  and  weaker 
than  the  0.16  ms  1  core  at  200m  observed  by  Hickey 
(1979).  During  this  period,  the  location  of  the  un¬ 
dercurrent  core  is  similar  between  NCOM-CCS  and  the 
CTW  model  (Fig.  5c),  which  shows  a  maximum  poleward 
velocity  of  0.06ms  V  Hie  strongest  mean  equatorward 
\  elocity  near  the  surface  in  the  I  lickev  ( 1 979)  observations 
is  0.08  m  s  1  over  the  outer  shelf,  which  is  much  weaker 
than  the  0.50  m  s  1  jet  during  the  sa  me  part  of  the  season 
in  NCOM-CCS.  As  discussed  further  in  section  5.  m- 
terannual  variability  is  a  likely  reason  for  the  difference 
between  the  1972  observations  and  2005  model  results. 

4.  Results 

Comparisons  with  observations  in  the  previous  sec¬ 
tion  indicate  that  the  CTW  and  NCOM-CCS  models  both 
capture  key  features  of  summer  circulation,  including 
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FlO.  4.  Comparison  of  CTW  modeled  alongshore  velocity  with  currenl-meler  measuremenis 
al  the  A1  n kk) ring  site  during  2004.  Model  lime  series  are  extracted  from  the  500* ni  isobath  at 
the  nominal  current -meter  depths 35.  100.  175.  and  5(X)m. 


equator  ward  How  at  midshclf  during  summer  and  \'ari- 
ability  in  the  California  Undercurrent  at  time  scales 
from  days  to  seasons.  I  he  weakest  aspect  of  the  models 
is  the  equatorward  shelf  break  jet.  which  is  absent  from 
the  (TW  model  but  present  over  the  Washington  slope 
in  NCOM-CCS  (Fig.  5).  To  gain  insight  into  the  linear 
dynamics' of  the  modeled  undercurrent,  the  CTW  model 
is  run  under  different  configurations  to  explore  the  sen¬ 
sitivity  of  the  2(X>5  time  series  and  the  seasonal  cycle  to 
friction,  number  of  modes,  local  wind  forcing,  and  the 
southern  boundary  condition  (section  4a).  Both  models 
are  then  used  to  determine  the  large-scale  structure  of 
the  APF  during  summer  (section  4b)  and  the  balance  of 
alongshore  momentum  in  the  northern  CCS  (section  4c). 

a .  CTW  model  sensitivity  analysis 

1)  Friction 

Previous  studies,  which  focused  on  the  shJf.  have 
shown  that  the  value  of  the  linear  frictional  coefficient 
can  i  nil  lienee  the  amplitude  and  phase  of  modeled  ve¬ 
locity  fluctuations  (Battisti  and  Hickey  NS4:  Chapman 
F)S7),  To  examine  how  friction  mod  ill  cs  the  response 
over  the  shelf  and  slope,  the  original  value  of  r  -  2.3  X 
10  4  ms  1  was  reduced  tor  —  1.0  X  10  4ms  *. and  also 
increased  by  the  same  amount  to  r  -  4.0  X  10  4ms  *. 
ITie  weaker  frictional  coefficient  enhances  the  modeled 
equatorward  velocity  oxer  the  midshelf,  which  reaches 

0.6  m  s  1  during  early  June  (Fig.  6a.  top),  much  larger 


than  the  observed  peak  velocity  of  0.3ms  1  (Fig.  2a). 
An  equivalent  increase  in  r  has  a  smaller  effect  on  the 
strength  of  equatorward  How  over  the  midshelf  (Fig.  6a. 
top).  At  300 m  oxer  the  slope,  decreasing  r  enhances 
poleward  Now  and  increasing/'  reduces  the  magnitude  of 
poleward  flow  during  late  summer  (Fig.  6b.  top),  with 
resulting  differences  <0.03  ms  *, 

2)  NUMBER  OF  MODES 

The  number  of  modes  used  x  anes  substantially  among 
studies  (Clarke  and  Van  Gorder  PFS6).  In  this  study, 
convergence  was  found  for  three  modes  and  more. 
During  mid-July -October  2005.  the  use  of  two  or  more 
modes  enhances  equatorward  floxv  oxer  the  shelf  (Fig. 
6a.  upper  middle)  and  enhances  polexvard  How  at  300  m 
over  the  slope  (Fig.  6b.  upper  middle ).  During  June  and 
early  July .  the  first  mode  dominates  the  response. 
Higher  modes  tend  to  be  locally  generated  because 
they  have  greater  x'alues  of  a„m  and  are  preferentially 
damped  as  they  traxel  away  from  the  forcing  region. 
Therefore  the  higher  modes  tend  to  be  more  important 
after  the  onset  of  strong  local  xvind  stress'  in  mid-July. 

3)  Local  wind  stress 

During  summer,  remote  wind  stress  contributes  to 
weather-band  variability  ox'er  the  Washington  shelf 
(Battisti  and  Hickey  1C)S4),  To  explore  the  role  of  re¬ 
mote  forcing  oxer  the  slope,  the  CTW  model  was  run 
with  local  wind  forcing  t*v  -  0  at  locations'  north  of42°N. 
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Fic».  5.  Sect  kins  of  alongshore  velocity  off  the  \VA  shelf,  (a)  Observed  geoslropltic  velocity  on  15  Sep  2(H15,  ref¬ 
erenced  lo  lUlJms”1  al  5(1(1  dhar  for  consistency  with  NC'OM-C’C'S  results.  Triangles  indicate  legations  of  (TD  casts, 
(h)  N('OM-C('S  results  for  15  Sep  2(105.  (c)  (T\V  model  results  for  15  Sep  2005.  (d)Tinie  average  of  NCOM-CCS 
results  for  23  Jul  26  Aug  2005.  (e)Tinie  average  of  CTW model  results  for  23  Jul  26  Aug  2(105.  Contour  intervals*  are 
().05ms  and  poleward  flow  is  shown  in  shades  of  gray  with  0.01  ms  1  intervals.  Labels  indicate  Ihe  location  and 
magnitude  of  the  strongest  cqualorward  velocity. 


'This  configuration  retains  the  relatively  strong  wind 
forcing  off  northern  California  and  removes  the  rela¬ 
tively  weak  forcing  north  of  Cape  Blanco,  Oregon  (Fig. 
Ic).  During  2005.  the  local  wind  stress  transitioned  to 
equatorward  over  one  month  later  than  usual  in  late 
May  and  remained  anomalously  weak  through  mid-July 
(Ilickcy  ct  al.  2006:  Kosro  et  al.  2006).  Over  the  shelf 
(Fig.  6a.  lower  middle),  remote  forcing  is  responsible  for 
the  strong  equatorward  How  (y  •*-  0.4  ms  l)  during  this 
early  period  of  weak  winds.  The  remotely  driven  shelf 
velocity  later  in  the  summer  during  August  is  weaker, 
with  more  frequent  poleward  How  events.  In  contrast  to 
the  shelf,  removing  the  local  wind  stress  has  little  effect 
on  the  flow  over  the  slope  at  300m,  aside  from  in¬ 
troducing  a  lag  of  several  days  (Fig.  6b.  lower  middle). 
The  difference  between  the  shelf  and  slope  can  he  un¬ 
derstood  by  considering  the  role  of  the  various  modes 


described  in  the  previous  section.  Removal  of  local 
equatorward  wind  stress  weakens  the  equatorward  flow 
associated  with  the  first  mode  over  both  the  shell*  and 
slope.  However,  the  suppression  of  local!)  generated 
higher  modes  results  in  weaker  equatorward  llow  over 
the  shelf  and  weaker  poleward  flow  over  the  slope. 
These  two  effects  counteract  each  other  over  the  slope 
but  reinforce  each  other  over  the  midshelf.  resulting  in 
greater  sensitivity  to  local  wind  stress  over  the  midshelf. 

4)  Southern  boundary  condition 

The  CTW  model  includes  signals  created  at  the 
southern  boundary  of  the  model  in  addition  to  those 
generated  by  wind  forcing  farther  north  along  the  west 
coast  of  North  America.  The  southern  boundary  con¬ 
dition  relies  on  the  assumption  that  fluctuations  in 
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CTW  mooel  velocity  sensitivity.  RN  location.  35m  (70m  isobath)  -  2005 


CTW  model  velocity  sensitivity  A1  location,  300m  (500m  isobath}  **  2005 


Fia  6.  Time  series  of  alongshore  velocity  showing  sensitivity  to  parameters  used  in  live  (TW 
model  at  (a)  35 m  over  lhc70-ni  isobath  at  the  RN  site  over  the  \VA  shelf  and  (b)  3(H)  ni  over  the 
500-m  isobath  at  the  A I  site  over  the  RC  slope.  In  each  row.  the  solid  black  is  the  base  run. 
Additional  runs  with  the  frictional  coefficients  r  -  1  X  10  4nis  1  (solid  gray)  and  r  -  4  X 
10  4nis  1  (dashed  gray)  (lop).  I  (dashed  gray)  and  2  modes  (solid  gray)  (upper  middle),  no 
wind  forcing  north  of425N  (solid  gray)  ( lower  middle  ),  and  no  energy  at  the  southern  boundary 
(solid  gray)  (bottom)  are  shown. 


adjusted  sea  level  at  San  Diego  are  associated  with  first- 
mode  coastal  trapped  waves.  I  o  show  the  effect  of  this 
southern  boundary  condition  on  the  model  results,  we 
usl  an  alternative  boundary  condition  of  f  |(().  /)  —  0. 
representing  zero  wave  energy  at  the  southern 
boundary. 

Over  the  midshclf,  the  southern  boundary  condition 
(Fig.  fia.  bottom)  docs  not  have  as  much  of  an  effect  on 
the  modeled  velocity  as  the  elimination  of  the  local  w ind 
stress  (Fig.  fia.  lower  middle).  However,  the  southern 
boundary  does  have  a  signilicant  effect  on  the  modeled 
velocity  over  the  slope  (Figs.  6b,  bottom).  Without  using 
the  San  Diego  sea  level  at  the  southern  boundary,  the 
3()()-m  \  cloeity  time  series  lacks  the  obse rved  persistence 
of  equatorward  flow  during  spring  and  poleward  (low 
during  late  summer  and  early  fall.  Poleward  flow  events 


still  occur  and  are  much  more  frequent  than  over  the 
shelf,  but  equatorward  flow  events  are  also  common 
during  late  summer. 

Seasonal-mean  sections  further  demonstrate  the  in¬ 
fluence  of  sea  level  variability  on  alongshore  velocity  in 
the  northern  CCS  (Fig.  7).  In  the  N COM -CCS  (Fig.  7. 
top)  and  the  CTW  model  with  San  Diego  sea  level 
variability  included  (Fig.  7.  middle),  poleward  flow  is 
a  maximum  of  — 0.0 5  m  s  1  at  300-400  m  duringsummer. 
intensities  near  the  coast  during  the  fall  and  winter,  and 
disappear;  during  spring  when  equatorward  flow  domi¬ 
nates.  Although  the  offshore  extension  and  vertical 
shear  of  equatorward  flow  are  more  realistic  in  NCOM- 
CCS  during  summer  (section  3;  Fig.  5).  the  development 
of  equatorward  flow  over  the  shelf  and  slope  during 
spring  and  the  intensification  of  poleward  flow  over  the 
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FIG.  7.  Seasonal  means  of  modeled  alongshore  velocity  at  47  N  during  (far  left)  January  March  2(H15,  (middle 
left)  April  June  20(15.  (middle  righl)  July  September  20(15.  and  (far  right)  October  December  20(15.  Results  are 
shown  for  the  (lop)  N COM- CCS  (middle)  CFVV  model,  and  (bottom!  CTW  model  with  no  wave  energy  at  the 
southern  boundary.  Contour  intervals  are  (1.05  ms  \  and  poleward  How  is  shown  in  shades  of  grav  with  (1.01  ms  1 
intervals. 


slope  during  summer  in  both  models  arc  consistent  with 
observational  studies  in  the  northern  CCS  (Hickey 
1070:  Collins  ct  al.  2003:  Thomson  and  Krassovski 
2010). 

In  the  case  with  no  CTW  energy  at  the  southern 
boundary  (Fig.  7.  bottom),  the  modeled  undercurrent  is 
weaker  in  the  summer  mean  (  0.02ms  l)  (Fig.  7).  In 
addition,  poleward  flow  over  the*  shelf  during  fall  and 
winter  as  well  as  equatorward  How  during  the  spring  are 
substantially  weakened  in  comparison  to  the  case  that 
includes  San  Diego  sea  level  variability.  Additional  runs 
that  include  winds  off  Baja  California,  south  of  San 
Diego,  showed  no  improvement  in  the  modeled  seasonal 
variations  in  the  northern  CCS  (not  shown).  Although 
CTW  dynamics  associated  with  wind  forcing  alone 
can  generate  weak  mean  poleward  How  over  the  northern 
CCS  slope  during  summer  (Fig.  7).  and  increase  polew  ard 


How  at  the  time  scale  of  weather  events  (Fig.  6b).  the 
magnitude  and  timing  of  its  seasonal  cycle  more  dosclv 
resemble  observations  and  NCOM-CCS  with  the  of 
San  Diego  sea  level  at  the  southern  boundary. 

/>.  Large-scale  structure  of  the  API •  during  summer 

The  structure  of  the  poleward  APF.  which  opposes 
the  prevailing  equatorward  wind  stress  during  summer, 
is  important  for  understanding  dynamics  associated  with 
the  CUC.  Although  the  sea  level  gradient  determines 
the  APF  near  the  surface,  the  APF  may  be  substantially 
modified  at  depth  bv  the  presence  of  alongshore  density 
gradients.  Cross-shelf  structure  may  also  be  present  due 
to  the  shell  topography  and  coastal  trapping.  In  this 
section,  we  use  the  NCOM-CCS  and  CTW  models  to 
gain  insight  into  the  structure  and  variability  ol  the  APF. 
Seasonal  differences  in  the  API  during  2(X>5  are  first 
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examined  near  the  coast,  where  tide  gauge  observations 
arc  available.  Analysis  of  the  modeled  API  is  then  ex¬ 
tended  offshore  to  the  continental  slope,  where  the 
models  provide  a  more  complete  view  of  its  timing  and 
spatial  structure . 

1  )  Y  ARI  ABILITY  IN  COASTAL  SEA  LEVEL 

Winter- to-su miner  differences  in  sea  level  are  com¬ 
pared  at  different  locations  along  the  coast  during 
January-March  and  15  June-15  September  2005  (Pig.  8a). 
The  latter  period  is  chosen  because  flow  transitions 
seasonally  from  equator-  to  poleward  at  300  m  over  the 
slope  (I  ig.  2c).  Differences  between  the  3-month  av¬ 
erages  are  examined  lust,  because  the  tide  gauge  ob¬ 
servations  cannot  be  used  to  estimate  the  absolute 
spatial  gradients  during  each  individual  season.  Be¬ 
tween  these  two  periods,  referred  to  as  “winter*'  and 
"summer"  conditions,  the  observed  sea  level  drops  to 
lower  values  along  the  entire  coastline  from  Point 
Conception  to  Vancouver  Island  (Pig.  8a),  consistent 
with  the  development  of  equatorward  gcostrophic  cur¬ 
rents  near  the  coast. 

Modeled  winter-to-summer  differences  in  sea  level 
follow  latitudinal  patterns  similar  to  those  of  the  ob¬ 
servations.  showing  larger  seasonal  decreases  in  the 
northern  CCS  (Fig.  8a).  Both  models  underestimate  the 
overall  magnitude  of  the  sea  level  decrease  at  all  lati¬ 
tudes.  Relative  differences  along  the  coast  in  NCOM- 
CCS  parallel  those  in  the  observations*  and  the  (T\V 
model  shows  a  weaker  alongshore  gradient  in  the  sea¬ 
sonal  differences.  To  reproduce  the  magnitude  of  the 
observed  sea  level  gradient  in  a  steady  barotropic 
model.  Hickey  and  Pola  (1983)  used  a  smaller  frictional 
coefficient  of/*  —  1  X  10  4ms  l.  Using  this  value  in  the 
CTW  model  results  in  a  larger  alongshore  gradient 
(Fig,  8a).  but  also  produces  unrealistically  strong 
equatorward  How  over  the  WA  shelf  (  up  to  0,6  m  s  l: 
Fig.  6.  top). 

Absolute  gradients  in  coastal  sea  level  from  the  nu¬ 
merical  modelscan  be  compared  with  obse nations  if  the 
tide  gauge  data  arc  added  to  the  mean  dynamic  height 
estimated  from  the  climatology.  Following  Hickey  and 
Pola  ( 1983),  mean  dynamic  height  data  at  the  coast  from 
Reid  and  Mantvhi  (1976),  referenced  to  5(X)dbar,  are 
interpolated  to  the  tide  gauge  locations.  During  winter 
2005,  the  N COM- CCS  and  CTW  models  both  show 
higher  average  sea  level  in  the  northern  CCS.  consistent 
with  the  observations  (I  ig.  8b).  During  summer,  coastal 
sea  level  in  NCOM-CCS  slopes  —  0. 10  m  downward 
from  PSL  in  the  south  to  Nil  farther  north,  consistent 
with  the  observed  alongshore  gradient.  Farther  north 
from  NI 1  to  LP.  the  obscrv cd  sea  level  gradient  indicates 
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coastal  sea  level  (2005  summer  -  winter  dillerence) 


PIG.  8.  Winter- lo-summer  differences  in  coastal  sea  level  from 
Point  Conception  to  Vancouver  Island  during  2(XJ5.  (a)  Difference 
of  seasonal  averages  from  adjusivd  lide  gauge  observations  (gray) 
and  N ( X)M- ( X N  ( Mac k  line).  CTW  (squares)  and  CTW  with  the 
drag  coefficient  reduced  by  a  factor  of  2 5  (triangles).  Winter  is  de¬ 
fined  as  January  March:  summer  is  defined  as  15  Jun  15  Sep. 
(b)  Summer  and  winter  averages  of  sea  level  for  the  three  models 
and  observations,  referenced  lo  the  latitude-averaged  mean  winter 
value.  Observ  ations  arc  estimated  by  adding  tide  gauge  data  to  the 
mean  dynamic  height  climatology  of  Reid  and  Manlyla  (1976). 

a  pok  ward  API  at  the  coast:  additional  data  in  northern 
CCS  indicate  that  this  pattern  is  robust  (  R.  McCabe  and 
B.  M.  Hickey  2013.  unpublished  manuscript).  In  con¬ 
trast  the  sea  level  gradient  in  NCOM-CCS  indicates 
a  weak  equatorward  APFat  the  coast  from  northern  OR 
to  southern  BC.  Smaller-scale  variability'  occurs  between 
LP  and  NB  in  NCOM-CCS  and  observations.  In  the 
CT  W  model,  coastal  sea  level  in  the  northern  CCS  slopes 
either  very  weakly  up-  or  downward  to  the  north,  de¬ 
pending  on  which  drag  coefficient  is  used.  Although  the 
absolute  gradient  in  the  coastal  sea  level  varies  between 
models  in  the  northern  CCS,  all  models  reproduce  the 
seasonal  development  of  the  CUC  (section  3).  As  will  lx 
shown  in  the  following  section,  this  is  due  to  the  im¬ 
portance  of  the  ba  roc  lime  component  of  the  APT  over 
the  slope. 

Modeled  intrascasonal  variability  in  the  APFovcr  the 
shelf  can  be  assessed  by  comparing  with  observed  vari¬ 
ability  in  alongshore  gradients  of  coastal  sea  level. 
Intrascasonal  fluctuations  in  the  observed  surface  API 
(  eh v.  where  h  is  sea  surface  height  at  the  coast)  be¬ 
tween  NH  and  NB  during  2005  are  more  correlated  with 
remote  alongshore  wind  stress  near  I  IB.  ov  er  400  km 
south  of  Newport,  than  with  local  wind  stress  (Fig,  9), 
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The  API7  is  directed  most  poleward  following  wind  re¬ 
laxations  or  reversals  to  the  south,  and  the  lags  of  thv 
maximum  correlation  with  wind  stress  along  the  coast 
suggest  a  response  that  travels  at  approximate!)  twice 
the  free-wave  (TW  phase  speed.  ITirs  propagation 
speed  of  the  pressure  response  is  consistent  with  ana¬ 
lytical  solutions  to  (7)  forced  by  periodic  wind  stress 
with  the  alongshore  structure  (Philander  and  Yoon 
F)S2).  Modeled  alongshore  gradients  in  coastal  sea  lex  el 
are  significantly  correlated  with  observations  in  the 
(TW  model  (CC  —  0.52:  Fig.  10a)  and  more  strongly 
correlated  with  observations  in  NCOM-CCS  (CC  - 
0.S7:  Fig.  10b).  Correlation  between  the  observed  and 
modeled  APF  at  the  coast  indicate  that  the  models  can 
be  used  to  understand  the  variability  of  the  APT  farther 
offshore  oxer  the  slope. 

2)  The  API  at  undercurrent  depths 

Analysis  of  the  modeled  APF  is  now  extended  from 
the  coast  to  the  slope,  where  the  CUC  ri  located.  To 
evaluate  the  slope  APF  at  the  same  alongshore  scale  as 
the  coastal  sea  level  gradients,  two  end  points  from  the 
same  latitudes  as  NB  and  Nil  are  selected.  Perturbation 
pressure  is  defined  as 

pL\.\\z.t)  p(e).  (M) 

where  JHz)  is  a  reference  pressure  calculated  from 
a  mean  density  profile,  which  does  not  contribute  to 
ho rizontal  grad ie nts. 

Shelf-to-slope  differences  in  the  modeled  APF  are 
significant  during  summer  (Fig.  MV).  Compared  with  the 
surface  APF  at  the  coast,  fluctuations  haxc  a  smaller 
magnitude  near  the  bottom  at  3IWm  ox  er  the  slope,  and 
there  is  no  significant  correlation  betxveen  the  API  at 
the  two  locations.  I  his  lack  of  correlation  is  likely  due  to 
higher-mode  wax  es.  The  spatial  and  temporal  x  ariabilitv 
of  bottom  pressure  at  3(Mm  during  this  period  shows 
several  poleward-propagating  signals  with  phase  speeds 
that  x  ary  betxveen  events  (Fig.  l()d).  An  example  during 
late  May  2005  is  highlighted  in  which  a  loxx  pressure 
signal  propagates  nor thxvard  at  approximately  twice  the 
phase  speed  of  a  free  first-mode  CTW.  followed  by 
a  high  pressure  signal  at  approximately  twice  the  phase 
speed  of  a  free  second-mode  CIAV.  During  this  event,  a 
poleward  APF  develops  oxer  the  slope  during  the  tran¬ 
sition  from  loxx  to  high  pressure  and  is  enhanced  bx  the 
poleward  decay  of  the  high  pressure  signal  (Figs.  lOc.d), 
Similar  patterns  can  be  found,  for  example,  in  the  time 
periods  leading  up  to  the  next  poleward  APF  event 
during  early  June  and  the  strongest  poleward  API  event 
during  mid-.luly. 


i(y)  vs.  Arj  (Newport  -  Neah  Bay) 


latitude 

Fig.  9,  (a)  Correlation  coefficient  between  observed  sea  level 
difference  (Nl I  M3)  and  CCMP  wind  stress  ai  various  latitudes. 
Filled  squares  indicate  .significance  at  the  95%  confidence  level, 
(b)  Lag  (days)  of  max  correlation. 


Flic  summer  means  of  modeled  perturbation  pressure 
and  SMI  are  used  to  further  examine  the  large-scale 
spatial  structure  of  the  APFxvith  depth  and  latitude  over 
the  slope  (Fig.  1 1).  South  of  42C>N.  the  mean  SSI  I  over 
the  slope  during  15  June— 15  September  is  consistent 
xvitli  a  large-scale  poleward  APF  at  the  surface  (Fig. 
1  la).  North  of  42°N  ox  er  the  slope,  higher  SSI  I  to  the 
north  indicates  a  mean  equatorxxard  APF  at  the'  surface. 
These  large-scale  patterns  in  the  surface  APF  arc  con¬ 
sistent  at  multiple  isobaths  oxer  the  slope.  Smaller-scale 
variability  (<l  latitude)  is  also  present  in  the  northern 
CCS.  Distributions  of  mean  bottom  pressure  oxer  the 
slope  (  Fig.  1  lb)  reveal  a  large-scale  poleward  APF  at  all 
latitudes  from  Point  Conception  to  Vancouver  Island. 
The  large-scale  API  at  the  bottom  is  weaker  at  isobaths 
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FIG.  IU.  Time  senes  of  observed  and  modeled  alongshore  pressure  gradients  duri ng  May 
September  21X15.  (a )  Observed  Nl  I  NB  sea  level  gradient  (gray)  and  (T\V  sea  level  gradient 
from  \VA  coast  segment  (black ).  (b )  Observed  NH  N B  sea  level  gradient  (gray]  and  sea  level 
gradient  calculated  from  the  same  coastal  locations  in  N ( X )M-( '( \S  ( black),  (c)  Sea  level  gradient 
calculated  from  the  same  coastal  locutions  in  NCOM-CCS.  with  bottom  pressure  gradient  at 
309  m  calculated  from  point  sat  the  same  latitudes  (dashed  gray).  The  bottom  pressure  gradient  at 
300 m  is  multiplied  bv  a  factor  of  5  for  clarity,  (d)  Bottom  perturbation  pressure  at  r  ~  3(19  m. 

from  Point  Conception  to  Vancouver  Island.  Dashed  lines  indicate  the  anticipated  propagation 
of  forced  Jirst-  and  second-mode  (TWs  (twice  the  free* wave  phase  speeds). 


deeper  and  farther  offshore.  The  horizontal  distribution 
of  perturbation  pressure  in  the  water  column  over  the 
3(l9-m  isobath  (Fig.  1  Ic)  shows  that  depths  >150m  arc 
characterized  by  lower  pressure  to  the  north  at  ah  lati¬ 
tudes.  indicative  of  a  poleward  APF  throughout  the  CCS 
at  depths  below  the  shelf  break. 

Hie  pressure  distribution  over  the  slope  indicates 
that  the  summer-mean  poleward  APFat  undercurrent 
depths  in  the  northern  CCS  is  associated  with  along¬ 
shore  density  gradients,  not  a  sea  surface  that  slopes 
downward  to  the  north.  Over  the  upper  slope,  maximum 
upward  displacement  of  the  s,  -  2(\.5  isopvcnal  occurs  at 
42  N  (Fig.  lid).  The  s,  -  2 5.5  isopvcnal  slopes  down¬ 
ward  to  the  north  at  42  N  and  lighter  isopyemris  (c.g.. 
s,  -  24.5)  also  slope  downward  to  the  north  in  the 
northern  CCS.  Therefore,  at  the  latitude  of  minimum 
sea  level  over  the  slope  and  latitudes  farther  north. 


alongshore  density  gradients  are  responsible  for  the 
poleward  API  at  undercurrent  depths.  Farther  south,  the 
alongshore  density  gradients  weaken  the  stronger  pole- 
ward  API  at  the  surface.  Similar  to  NCOM-CCS.  the 
maximum  density  anomaly  in  the  CTW  model  occurs 
near  42  N  (Fig.  11c).  While  the  magnitude  of  these 
anomalies  are  about  an  order  of  magnitude  greater  in 
NCOM-CCS.  the  pattern  of  upward  isopvcnal  dis¬ 
placement  centered  at  42  N  is  consistent  with  CTW 
dy  namics.  This  upward  displacement,  most  pronounced 
in  the  s,  -  26.5  isopvcnal.  docs  not  appear  to  result  from 
offshore  forcing  because  it  is  much  weaker  over  the 
1055- m  isobath  in  NCOM-CCS  (Fig.  1  If). 

The  CTW  model  can  be  used  to  evaluate  contribu¬ 
tions  to  the  API  from  alongshore  wind  forcing  along  the 
U.S.  West  Coast  and  remote  signals  farther  south.  The 
cross-shore  structure  of  the  large-scale  APF  is  clearly 
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three  isobaths,  (b  )  NCOM-CCS  lx>ttoni  pert urKition  pressure  over  three  isobaths.  normalized  to  units  of  dynamic 
height,  (c)  NCOM-CCS  perturbation  pressure  over  the  309-ni  isobath,  (d)  NCOM-CC'S  velocity  y  and  s,  (contours) 
over  the  309-ni  isobath.  (e)  CTW  alongshore  y  (colors)  and  density  anomaly  (contours)  over  the  3(W-ni  isobath. 
(I)  NCOM-CCS  velocity  y  and  s,  (contours)  over  the  1055- m  isobath.  Note  the  smaller  contour  interval  for  CTW 
density  anomaly  compared  with  NC'OM  C'C'S  Sf. 


defined  in  the  CTW  model,  because  the  bathymetry  is 
uniform  in  the  alongshore  direction  within  each  seg¬ 
ment.  During  15  June-15  September  2005.  the  poleward 
API  in  the  northernmost  segment  of  the  CTW  model  is 
strongest  near  the  bottom  at  the  shelf  break,  and  an 
equatorward  API  is  present  over  the  shelf  near  the  coast 
(Pig.  I2a).  Farther  south  in  the  CTW  model,  at  the  17  N 
Washington  segment,  the  strongest  poleward  API  oc¬ 
curs  closer  to  shore  over  the  shelf  (Pig.  1 2d ) .  At  the 
40. S  N  I  IB  segment,  closer  to  the  strongest  wind  forcing, 
the  poleward  API7  is  strongest  at  the  surface  and  near 
the  coast  (Fig.  I2g). 

The  stronger  poleward  APT  at  the  depth  over  the 
slope  (Fig.  12a)  again  demonstrates  the  importance  of 
alongshore  density  gradients  at  locations  north  of  the 
strongest  wind  stress  in  the  CCS.  Unlike  the  NCOM- 
CCS  summer  mean  (Pig.  I  la),  the  CTW  model  does  not 
have  an  equatorward  surface  APT  over  the  northern 
CCS  slope,  which  may  be  associated  with  the  equator- 
ward  shelf  break  jet  that  is  absent  from  the  Cl AV  model. 
However,  the  NCOM-CCS  and  CTW  models  both  in¬ 
dicate  enhancement  of  the  poleward  API7  at  depth  over 
the  slope  by  alongshore  density  gradients.  In  the  fol¬ 
lowing  section,  it  will  be  shown  that  the  modeled  pole- 
ward  APT  over  the  northern  CCS  slope  is  largely 


balanced  by  a  Coriolis  force  associated  with  onshore 
flow  at  time  scales  from  weeks  to  months. 

<\  Alongshore  momentum  balance  over  the  northern 
CCS  slope 

In  a  linear  model,  the  alongshore  momentum  balance 
away  from  the  frictional  boundary  layers  is  described  by 

y,  +  C/y  (l(» 

The  CTW  model  indicates  that  during  summer  the  mean 
APT  r0  Vv  is  primarily  balanced  by  the  mean  Coriolis 
force  fit  (Pigs.  I2h.c.h),  which  is  much  greater  than  the 
mean  acceleration  y,  (Figs.  I2c.f.i).  In  this  section,  tem¬ 
poral  variability  in  (10)  is  examined.  The  validity  of 
a  linear  alongshore  momentum  balance  is  assessed,  and 
it  is  shown  that  the  dominant  terms  in  the  alongshore 
momentum  balance  are  dependent  on  the  time  scale  of 
interest,  with  the  Coriolis  force  largely  balancing  the 
API7 at  >20-dav  timescales. 

The  validity  of  the  linear  momentum  balance  over  the 
slope  can  he  examined  using  the  nonlinear  NCOM-CCS. 
Io  remove  mcsoscale  variability  and  focus  on  the 
alongshore  scale  of  the  northern  CCS  slope.  NCOM- 
CCS  time  series  have  been  averaged  over  44.5  — 48.5  N 
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FIG.  12.  Seasonally  averaged  momentum  balance  terms  from  the  (TVV  model  during  15  Jun-15  Sep  2005:  (a)  mean 
APF.  (b)  mean  Coriolis  force,  and  (c)  mean  acceleration  at  the  ‘IT  segment  off  the  British  Columbia  coast  (see  Fig.  I ). 
(d)  (f)  As  in  (a)  (c).  but  for  the  WA  segment  (see  Fig.  I),  (g)  (i)  As  in  (a)  (c).  but  for  the  I  IB  segment  off  the 
northern  coast  of  California.  Contour  intervals  are  I  x  10  ?  iti  s  2  with  positive  values  shaded  at  0.2  x  10  1  ms  2 
intervals. 


and  3  grid  points  offshore  of  the  309-m  isobath.  To 
separate  variability  at  time  scales  of  several  weeks  or 
longer  from  weather-band  variability,  low-frequency 
time  series  have  been  created  by  filtering  with  a  20-dav 
1  fanning  window. 

In  the  low-frequency,  spatially  averaged  NCOM- 
CCS  time  series,  the  undercurrent  accelerates  and 
shoals  during  summer  (Fig.  1 3a).  The  APF  and  Coriolis 
force  resemble  each  other  closely  in  timing,  magnitude, 
and  depth  structure  (Figs.  13b.c).  Although  there  is 
a  mean  haroclinic  depth  structure  (Fig.  lie),  a  surface- 
intensified  poleward  APF  is  present  during  July  (Fig. 
13b).  Acceleration  is  generally  an  order  of  magnitude 
smaller  than  the  APF  and  Coriolis  force  at  these  longer 
time  scales  (Fig.  13d).  Agreement  between  acceleration 


and  the  residual  of  the  Coriolis  force  and  APF  away  from 
the  surface  boundary  layer  indicates  that  the  linear  bal¬ 
ance  is  valid  over  the  slope,  especially  at  100-300  m  (Figs. 
I3d.c). 

At  a  typical  undercurrent  depth.  50  m  above  the  bot¬ 
tom  boundary  layer  at  259  m.  the  dominant  tc  mis  of  the 
momentum  balance  depend  on  the  time  scale  analyzed. 
At  shorter  time  scales  (<20-dav  periods),  the  APF  is 
generally  smaller  than  the  acceleration  and  Coriolis 
force  terms  (Fig.  14).  The  timing  and  magnitude  of  the 
momentum  balance  terms  arc  consistent  between  the 
N  COM  -CCS  and  CTW  models.  Correlations  between 
models  arc  0.55.  0.50.  and  0.49  for  the  acceleration. 
Coriolis,  and  APF  terms,  respect  holy,  with  a  1-dav  lag  in 
the  APF  term. 
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Fici.  13.  Velocity  and  linear  momentum  balance  terms  (20-day  period  and  greater)  over  the 
northern  (X'S  slope  (44.5*  4S.5”M  during  2005  from  NCOM*(X'S:  (a)  alongshore  velocity, 
(h)  APF.  (c)  Coriolis  force,  (d)  alongshore  acceleration,  and  (e)  residual  of  the  APF  and 
Coriolis  force.  Note  the  color  scale  difference  between  (b)  (c)  (XIO  h)  and  (d)  (e)  ( x|()  7). 


At  longer  time  denies,  the  Coriolis  force  largely  bal¬ 
ances  the  APF  and  leads  the  API*  by  2  days  in  NCOM- 
CCS  (Fig.  1 5a).  In  the  (TW  model,  the  Coriolis  force 
and  APF  arc  weaker  and  not  as  strongly  correlated  (I  ig. 
15b).  I  lowever.  the  seasonal  timing  of  the  APF  is  similar 
in  the  two  models,  with  maxima  in  the  poleward  APF 
during  late  April  and  mid-Julv.  In  both  models,  the 
poleward  APF  is  persistent  during  the  period  June- 
August.  The  acceleration  term  (not  shown)  is  weak  and 
lacks  a  persistent  summer  mean  (see  Figs.  12,  13).  With 
wind  forcing  only  [f  i((). /)  —  (>|.  the  Coriolis  term  and  the 
APF  are  weaker,  but  onshore  flow  and  the  poleward 
APF  are  still  persistent  during  summer  (Fig.  15c).  F\cn 
though  the  mean  poleward  velocity  during  summer  is 
weaker  in  the  wind-only  case  (Figs.  6b.  7).  the  seasonal 
development  of  the  poleward  APF  over  the  northern 
CCS  slope  is  a  robust  part  of  each  model  run  and 
therefore  is  largely  associated  with  wind  forcing  along 
the  U.S.  West  Coast. 

5.  Summary  and  discussion 

Hits  study  uses  two  numerical  models  with  different 
levels  of  complexity  to  examine  the  dynamics  of  the 
poleward  CUC  and  APT .  Time  series  observations  of 
alongshore  velocity  over  the  continental  slope  allow  for 
a  detailed  assessment  of  modeled  variability  over  time 
scales  from  days  to  seasons.  Processes  influencing  the 


CUC  and  APF  are  illustrated  schematically  in  Fig.  16, 
and  major  results  arc  summarized  below: 

*  Seasonal  development  of  a  persistent  CUC  duringlatc 
summer  in  the  northern  CCS  is  primarily  dependent 
on  incoming  energy  from  sea  level  variability  applied 
south  of  the  region  of  strong  coastal  wind  forcing, 
while  event-scale  fluctuations  over  the  northern  CCS 
slope  at  CUC  depths  are  primarily  forced  by  remote 
coastal  wind  events  (section  4a). 

«  The  APF  over  the  continental  slope  has  barotropic 
and  baroclinic  components.  At  time  scales  from  weeks 
to  months  in  the  northern  CCS,  the  baroclinic  API 
provides  a  poleward  force  at  CUC  depths,  which  is 
balanced  by  the  Coriolis  force  associated  with  onshore 
How .  In  the  central  and  southe rn  parts  of  the  CCS,  the 
baroclinic  component  weakens  the  barotropic  APF  at 
CUC  depths  (sections  4b-c ). 

An  unanticipated  result  of  this  study  is  the  de¬ 
pendence  of  a  persistent  late-summcr  CUC  on  the  use 
of  the  San  Diego  sea  level  at  the  southern  boundary. 
Low  adjusted  sea  level  at  San  Diego  during  spring  is 
associated  with  strong  equatorward  flow,  while  high 
adjusted  sea  lc\el  is  associated  with  a  poleward-llowing 
inshore  countercurrent  (Reid  and  Ylantyla  1 l)76).  This 
near-surface  poleward  flow  is  evident  in  the  mean 
NCOM-CCS  sea  surface  height  field  during  summer 
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Fl(l.  14.  Modeled  linear  moment  uni  balance  terms  (<20-day  periods  only)  during  2005  from 
the  NCOM-CCS  and  CTW  models:  (a)  acceleration,  (b)  Coriolis  force,  and  (c)  APF.  Time 
series  are  for  50 m  from  the  bottom  over  the  309-m  isobath.  NCOM-CCS  lime  series  arc  av¬ 
eraged  over  44.5:  ~4S.5  N.  and  CTW  lime  series  arc  from  the  WA  coast  segment. 


(Tig.  la).  Model  results  have  shown  that  remote  forcing 
from  the  equatorial  Pacific  contributes  to  the  annual 
cycle  of  sea  level  in  the  CCS  (Pa res -Sierra  and  OTJrien 
19X9).  More  recently.  Shulman  et  al.  (21)10)  linked  ob¬ 
servations  of  anomalously  high  coastal  sLa  level  during 
summer  2006  to  propagating  signals  originating  near  the 
equator  in  the  global  NCOM  model  1  he  importance  of 
the  southern  boundary  at  seasonal  time  scales  suggests 
a  possible  dynamic  link  between  poleward  Hows  often 
considered  to  be  separate:  the  undercurrent  as  far  north 
as  British  Columbia  and  the  inshore  countercurrent  in 
the  southern  CCS. 

Ill  is  result  contrasts  with  the  study  of  Chapman 
( 19X7).  who  found  that  including  San  Diego  sea  level  as 
a  southern  boundary  condition  weakened  the  correla¬ 
tion  with  observed  bottom  pressure  over  the  northern 
CA  shelf.  In  the  present  study,  the  southern  boundary* of 
the  CTW  model  primarily  influences  the  northern  CCS 
at  a  time  scale  of  months,  which  (Tiapman  ( 1 9X7)  did  not 
analyze.  It  is  possible  that  smaller-amplitude  weather- 
band  fluctuations  at  San  Diego  are  generated  more  lo¬ 
cally  than  the  lower- frequency  fluctuations  and  are 
therefore  associated  more  with  higher-mode  CTWs  that 
do  not  propagate  far  north.  It  is  also  possible  that 
the  shorter-period  weather-band  fluctuations  are  more 
susceptible  to  scattering  by  alongshore  variations  in 


topography,  which  most  likely  occurs  whe  re  the  coastline 
bends  sharply  near  Point  Conception  (south  of  PSL:  Fig. 
la).  Johnson  ( 1991 )  shows  that  scattering  can  be  limited 
cither  if  stratification  is  strong  or  if  frequencies  are 
much  less  than  AVa.  where  a  is  the  bottom  slope.  For 
characteristic  values  near  Point  Conception.  A;2  - 
10  4s  1  and  a  -  5  x  10  3:  this  corresponds  to  a  period 
of  —1.5  days,  which  is  comparable  to  the  time  scales  in 
the  weather  band. 

Outside  of  the  northern  CCS.  undercurrent  dynamics 
may  be  more  strongly  influenced  by  processes  such  as 
wind  stress  curl  that  arc  not  present  in  the  CTW  model 
dynamics  or  implicitly  included  in  the  southern  bound¬ 
ary  condition,  McCreary  et  al.  (19X7)  and  Battecn  ( 1997) 
found  that  wind  stress  curl  creates  a  stronger  and  shal¬ 
lower  poleward  undercurrent.  However,  remote  wind 
stress  curl  is  not  thought  to  influence  coastal  currents 
in  the  same  way  as  remote  alongshore  wind  stress 
(McCreary  et  al.  19X7).  Poleward  flow  induced  by  the 
alongshore  gradient  of  wind  stress  curl  weakening  to  the 
south  may  influence  remote  locations  (Ocv  1999).  but  it 
is  unclear  how  this  mechanism  would  generate  poleward 
flow  in  the  northern  CCS.  Wind  stress  curl  is  relatively 
weak  in  the  northern  CCS  and  strongest  off  northern 
California  and  in  the  Southern  California  Bight  (Bakun 
and  Nelson  1991). 
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FKi.  15.  Modeled  APF  and  Coriolis  force  terms  (20-day  period  and  greater)  over  the 
northern  CCS  slope  during  2(X)5.  (a)  NCOM*C('S.  (b)  (TW  model,  and  (c)  (TW  model  run 
with  no  wave  energy  at  the  southern  boundary.  All  lime  series  are  from  259-ni  depih  over  the 
309- ni  isobath.  as  in  Fig.  14.  NCOM-CC'S  time  series  arc  averaged  over44.5*-48.59N.  and  (TW 
time  series  are  from  ihe  \VA  coast  segment. 


Another  factor  missing  from  the  (TW  model  is  the  b 
effect,  which  in  idealized  studies  has  been  found  to 
strengthen  the  poleward  undercurrent  and  APF  at  sea¬ 
sonal  time  scales  (Philander  and  Yoon  1982:  Suginoharu 


and  Kitamura  1984).  Observations  in  the  southern  CCS 
have  revealed  offshore  propagation  of  poleward  How 
structures  consistent  with  baroclinic  Rossbv  waves 
(  Todd  ct  al.  2011).  However,  it  remains  unclear  how 
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Fig.  16.  Schematic  of  ]  Processes  influencing  the  ClJCand  APF  during  summer.  The  term  hv 
represents  the  alongshore  pressure  gradient  force  associated  with  sea  level  hCUMt  represents  sea 
level  at  the'  coast.  rv  represents  the  alongshore  density  gradient  over  the  slope,  and  t  'represents 
Ihe  alongshore  component  of  wind  stress.  Solid  red  arrows  represent  currents.  Mack  dashed 
arrows  represent  the  APF.  and  thick  blue'  lines  represent  isopvenals  intersecting  a  vertical  plane 
over  the  slope. 
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important  the  b  effect  is  at  seasonal  and  shorter  time 
scales  farther  north  of  Washington  and  British  Colum¬ 
bia.  At  the  annual  time  scale,  Rossbv  wave  propagation 
only  occurs  south  of  a  critical  latitude,  estimated  to  be 

42  N  in  the  CCS  (Clarke  and  Shi  1991). 

I’opographic  rectification  is  another  process  that  is 
absent  from  the  CTW  model  but  could  he  present  in 
NCOM-CCS.  In  the  presence  of  stratification.  the  rec¬ 
tification  of  oscillating  How  over  topography  that  varies 
m  the  along-shelf  direction  produces  a  mean  poleward 
How  concentrated  over  the  upper  slope.  The  modeling 
study  of  Brink  (2011)  shows  that  parameters  represen¬ 
tative  of  the  northern  CCS  result  in  a  0.1  ms  1  mean 
poleward  flow.  In  addition,  topographic  rectification  can 
also  lead  to  a  large-scale  gradient  in  coastal  sea  level  that 
slopes  downward  in  the  direction  of  wave  propagation 
(Brink  2010).  It  is  therefore  possible  that  topographic 
rectification  could  explain  the  stronger  poleward  flow 
and  coastal  sea  level  gradient  in  NCOM-CCS  compared 
with  the  (TW  model  during  summer,  although  the 
—40km  topographic  length  scale  is  only  coarsely  re- 
solvod  by  the  9-km-rcsolution  NCOM-CCS.  Because 
the  topographic  rectification  theory  provides  scaling  for 
the  mean  flow  ,  it  is  not  immediately  clear  w  hether  it  can 
explain  the  temporal  variability  that  lias  been  the  pri¬ 
mary  focus  of  the  present  study.  Seasonally,  the  sub¬ 
surface  poleward  How  associated  with  topographic 
rectification  is  expected  be  strongest  when  variability  in 
the  local  wind  stress  is  highest.  In  the  northern  CCS. 
the  standard  deviation  of  daily  mean  wind  stress  is 
strongest  during  December  and  reaches  a  minimum 
during  August  (Hickey  1979).  Therefore,  although  to¬ 
pographic  rectification  cannot  be  ruled  out  as  a  potential 
mechanism,  it  docs  not  appear  to  explain  the  intensi¬ 
fication  of  poleward  flow  over  the  upper  slope  from 
spring  to  late  sum  me  r. 

The  strength  of  the  CUC  also  varies  on  interannual 
time  scales,  which  may  explain  differences  between 
seasonal  averages  of  the  NCOM-CCS  results  from  2(X)3 
and  the  Hickey  (  1479)  current-meter  observations  from 
1972  (section  3).  Summer  1 972  coincides  with  a  strong  HI 
Nino,  although  it  is  not  associated  with  significant  warm 
temperature  anomalies  in  the  northern  CCS  like  other 
I  I  Nino  years  (Smith  et  al.  2001)  However,  current- 
meter  measurements  indicate  much  stronger  poleward 
velocity  oxer  the  outer  shelf  and  slope  during  late 
summer  1972  than  late  summer  1981  (Hickey  1989), 
Mean  late-summcr  poleward  flow  at  200m  over  the 
upper  slope  is  —0.03  ms  1  during  1981.  compared  with 

0.13  ms  1  during  1972.  Mean  late-summcr  current- 
meter  observations  from  1981  show  near-surface 
cquatorward  flow  of  about  0.20ms  1  (Ilickev  1989). 
compared  with  0.08  m  s  1  during  the  1972  observations. 


Over  midshclf,  at  the  80-m  isobath,  late-summer  ob¬ 
servations  from  both  1972  and  1981  indicate  that  pole- 
ward  flow  extends  throughout  nearly  the  entire  water 
column.  1  Iowever.  neitherobserv ationsnor  model  results 
show  mean  poleward  flow  in  the  midshelf  water  column 
during  late  summer  2003  (section  3).  The  presence  of 
stronger  cquatorward  flow  over  midshelf  during  2(X)5 
than  other  years  may  be  related  to  anomalous  local  wind 
stress,  which  transitioned  cquatorward  about  l  month 
later  than  usual  off  the  Washington  coast  but  was  stronger 
than  usual  during  the  later  pan  of  the  summer  (Ilickev 
et  al.  2(XKk  Kosro  et  al  2(X)6).  Sensitivity  of  alongshore 
currents  to  local  wind  stress  in  the  (TW  model  is  higher 
over  the  shelf  than  the  slope  (section  4a).  It  is  therefore 
likely  that  in  ter  annual  variability  of  the  ClTC  in  the 
northern  CCS  is  more  strongly  influenced  by  remote 
forcing  than  local  wind  stress  and  also  influenced  by 
b-planc  dynamics  as  discussed  above. 

Hie  baroclinic  structure  of  the  API  presented  in 
this  study  is  consistent  with  the  idealized  model  of 
McCreary’  (1981),  in  which  the  undercurrent  deepens  at 
higher  latitudes  and  reaches  a  maximum  amplitude  at  the 
latitude  where  isopycnal  surfaces  arc  shallowest  (Fig. 
t  Id).  The  presence  of  a  poleward  API'  at  depth  over  the 
slope  does  not  depend  on  the  southern  boundary*  condi¬ 
tion  of  tliL  CTW  model  (Figs.  1 3b-c)  or  a  strong  polcwa rd 
API  near  the  coast  (Figs.  I2a-c).  McCreary  and  (Turn 
( 1983)  show  that  this  structure  of  the  APT  is  robust  ev  on 
when  the  inclusion  of  a  continental  shelf  weakens  the 
undercurrent.  The  same  basic  structure  is  shown  in  the 
two-layer  model  of  Yoshida  { 1967)  in  which  the  sea  level 
reaches  a  minimum  at  the  northern  edge  of  the  wind 
stress  forcing  and  the  shallowest  depth  of  the  thermocline 
occurs  at  the  same  latitude.  In  the  present  study*  over  the 
CA  slope,  the  API  is  poleward  through  the  entire  water 
column  and  alongshore  density  gradients  weaken  the 
poleward  API  with  depth  (Figs,  lla-c),  consistent  with 
previous  studies  that  focused  on  these  latitudes  (Mellor 
1986;  Pringle  and  Dover  2(X)9). 

Model  hindcast  results  in  the  present  study  indicate 
that  the  largest  vertical  isopycnal  displacements  over  the 
slope  occur  north  of  the  strongest  cquatorward  wind 
forcing  (Fig.  II),  in  agreement  with  past  idealized 
models  (Yoshida  t%7:  McCreary  1981).  Shipboard  ob¬ 
servations  also  suggest  a  maximum  vertical  isopycnal 
displacement  over  the  slope  that  is  offset  to  the  north  of 
the  strongest  upwelling-favorablc  winds  (Muraki  1974). 
Vertical  motions  of  isopycnal s  over  the  slope  are  asso¬ 
ciated  with  higher  (TW  modes,  which  are  damped  more 
quickly  than  the  lirst-niode  waves  that  dominate  the 
coastal  sea  level  response.  Higher  background  strati¬ 
fication  to  the  north  may  also  contribute  to  the  baro- 
ctinic  API  by  restricting  the  upward  displacement  of 
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isopycnals.  In  the  (TW  model,  onshore  flow  associ¬ 
ated  with  the  poleward  API’  at  the  shell  break  trans¬ 
fers  mass  from  the  poleward  undercurrent  to  the 
cquatorward  shelf  jet  in  the  northern  CCS.  In  NCOM- 
CCS.  the  baroelinie  structure  of  the  APF  and  associated 
cross-shore  flow  near  the  surface  are  complicated  by 
the  presence  of  equator  ward  flow  farther  offshore  o\cr 
the  shelf  break  (Fig.  5).  Observations  have  shown  that 
the  cquatorward  jet  migrates  offshore  throughout  the 
summer  (Hickey  1981)),  which  may  influence  the  mag¬ 
nitude  and  structure  of  the  poleward  APF  near  the 
surface  in  the  northern  CCS.  The  contribution  of  density 
gradients  to  the  poleward  APF  over  the  slope,  demon¬ 
strated  here  in  the  northern  CCS.  is  expected  to  he  dy¬ 
namically  important  in  the  poleward  latitudes  of  other 
eastern  boundary  upwclling  systems. 

Tliis  study  presents  a  comprehensive  analysis  of  sea¬ 
sonal  a  nd  e\  e  nt-scale  flow  ove  r  the  continental  slope  of  the 
CCS.  Models  are  tested  against u  uniquclong-tc  mi  dataset 
collected  over  the  northern  CCS  slope.  Our  results  high¬ 
light  the  importance  of  remote  forcing  in  all  parts  of  the 
COS.  which  provides  motivation  for  treating  such  systems 
asa  unified  whole  in  future  modeling  studies.  The  presence 
of  strong  variability  at  multiple  time  scales  demonstrates 
that  the  collection  of  in  situ  data  over  the  continental  slope 
should  be  a  critical  component  of  future  dynamical  studies 
of  poleward  undercurrents  in  upw  clling  regions. 
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